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Results of shell-model calculations for lower fp-shell nuclei show that SU{3) symmetry breaking 
in this region is driven by the single-particle spin-orbit splitting. However, even though states of 
the yrast band exhibit 517(3) symmetry breaking, the results also show that the yrast band B{E2) 
values are insensitive to this fragmentation of the SU{3) symmetry; specifically, the quadrupole 
^^ ' collectivity as measured by B{E2) transition strengths between low lying members of the yrast 

^"—^ , band remain high even though SU{3) appears to be broken. Results for 44,46,48j.^ ^^^ '^^Cr using 

the Kuo-Brown-3 two-body interaction are given to illustrate these observations. 
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PACS number (s) : 21. 60. Cs, 21.60.Fw, 27.40.-l-z 



I. INTRODUCTION 



SU{2i) is a special algebraic structure because it is the compact symmetry group of the three-dimensional isotropic 
harmonic oscillator [y which is a good first-order approximation to any attractive potential. This applies in nuclear 
^^' physics and is the underpinning to the Elliott S't/(3) model B]. In the latter case the highest symmetry group is 
f^ , SU{Kfl), where J7 denotes the degeneracy of the spatial degrees of freedom and k counts the number of internal 
f^ • degrees of freedom (for example, k — 2 for identical spin i particles and k = 4 for a spin-isospin system). SU{i) 
enters in this picture through a reduction of SU{Kil) into spatial [S'C/(0)] and spin or spin-isospin degrees of freedom 
[SU{k)], namely, SU{kV,) D SU{Q)(>^ SU{k), followed by a reduction of the spatial degrees of freedom through SU{3) 
to 5*0(3); that is, SU{n) D SU{3) D 50(3). Interactions that are not functions of the 5C/(3) generators induce 
O , 5t/(3) symmetry breaking. The spin-orbit interaction, which is needed for a correct description of shell and subshell 
closures |^, is an example of a one-body SU{3) symmetry breaking interaction while the pairing interaction, which 
is required for a correct description of binding energies y, is an example of a two-body SU{3) symmetry breaking 
interaction. 

'k^ ' It is well known that SU{3) is a very useful symmetry in the lower sd-shell [||. This is most easily understood 

'Oj . by noting that the leading irreducible representation (irrep) of SU{3) normally suffices to achieve a good description 
d of the low-lying eigenstates of these nuclei. In the lower /p-shell, however, leading SU{3) irreps do not provide 
satisfactory results for low-lying eigenstates. Beyond the /p-shell, the concept of pseudospin symmetry B allows one 
to identify another so-called pseudo SU{3) structure that again yields a good description of low-lying eigenstates of 
strongly deformed nuclei . Questions that remain regarding the lower /p-shell are: What parts of the interaction 
are responsible for the 5J7(3) symmetry breaking? Is it the one-body part, the two-body part, or a combination of 
these two? And if it is a combination, to what extent does each interaction contribute to SU(3) symmetry breaking? 
Also, what is the effect of the SU{3) symmetry breaking on the electromagnetic transition rates? Enhanced B{E2) 
transition rates ||7| are normally considered to be a good indicator of quadrupole collectivity and the SU{3) structure of 
the corresponding initial and final states. It has been suggested that strong B{E2) values may survive an "adiabatic" 
mixing of SU{3) irreps due to quasi-5C/(3) dynamical symmetry [||. A signature for this type of mixing is B{E2) 
values that are similar to those obtained when the SU{3) symmetry is good. Is the SU{3) symmetry breaking in the 
lower /p-shell adiabatic? 

In this paper we show for lower /p-shell nuclei that whereas the spin-orbit interaction is the primary driver of 
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SU{3) symmetry breaking the B{E2) values between the first few yrast states remain strong, signahng an adiabatic 
mixing of SU{3) irreps. The reahstic monopole-corrected Kuo-Brown-3 two-body interaction [^ is used in calculations 
for 44,46,482ij g^jjj "^^Cr with single-particle energies corresponding to realistic spin-orbit splitting. The spectrum of 
the second-order Casimir operator C2 of SU{3) is used as a measure for gauging the SU{3) fragmentation along the 
yrast band of these nuclei. The results show that the spin-orbit splitting is the primary cause for SU{3) symmetry 
breaking; the leading SU{S) irrep regains its importance as the spin-orbit splitting is turned off. A similar recovery 
of the SU{3) symmetry has been reported in the case of '^'^Ti with degenerate /i — pi shells [ |lO[ . 

To fix the notation, in the following section a parameterization of the Hamiltonian in terms of one-body spin-orbit 
and orbit-orbit single-particle interactions, as well as a general two-body interactions, is given. In our applications of 
the theory, the realistic Kuo-Brown-3 interaction is chosen for the two-body interaction f^. Computational methods 
used in the analyses are discussed in the third section. This is followed by characteristic results for '^'^Ti, ^^Ti, '^^Ti, 
and ^^Cr in the fourth section. A conclusion that recaps outcomes is given in the fifth and final section. 

II. INTERACTION HAMILTONIAN 

The one- plus two-body Hamiltonian is used in standard second-quantized form: 



H = J2'^^a+a, + jJ2 ^ki.tja+a^j 



4 ^^ ■r..,.J-~^ ^^jakai. 

i i.j,k,l 

The summation indexes range over the single-particle levels included in the model space. We only consider levels 
of the /p-shell which have the following radial (n), orbital {I) and total angular momentum (j) quantum numbers: 

nlj = < 0/7 , 0/5 , Ipi , Ipj. \ . In what follows the radial quantum number (n) is dropped since the Ij labels provide a 

unique labelling scheme for single-shell applications. It is common practice to replace the four single-particle energies 
Ei in the /p-shell by the P and I ■ s interactions: ^^ eiO^ai — > e{ni — aiU ■ Si — (3ilf): where e is the average binding 
energy per valence particle, rii counts the total number of valence particles, and a and (3 are dimensionless parameters 
giving the interaction strength of the P and l-s terms. For realistic single-particle energies used in the KB3 interaction 
(|l|), these parameters are e = 2.6 MeV, (3 = 0.0096, Up = 1.3333, and a/ — 1.7143. The small value of j3 signals small 
f splitting (|). 

A significant part of the two-body interaction, Vki^ij , maps onto the quadrupole-quadrupole, Q ■ Q, and the pairing, 
P, interactions. Since Q ■ Q can be written in terms of SU{3) generators, it induces no SU(3) breaking and hence 
serves to re-enforce the importance of the Elliott model g , when the pairing interaction mixes different SU{3) irreps. 
In this analysis the two-body part of the Hamiltonian, Vki^ij, is fixed by the Kuo-Brown-3 (KB3) interaction matrix 
elements while the single-particle energies, Si, are changed as described below. 

The following single-particle energies are normally used with the KB3 interaction Q: 

KB3 [MeV]:epi=4, ep., = 2, (1) 

For the purposes of the current study, it is important to know the single-particle centroids of the p and / shells. 
For example, the energy centroid of the p shell is given by: 

epj_ dim(pi ) + Sp^ dim(p3) 

_ 2 ^ 2 ^__ 

dim(pj.) + dim(p3) 

In what follows, we label by KB3p-f that Hamiltonian which uses the KB3 two-body interaction with single-particle 
p- and /-shell energies set to their centroid values: 

KB3p_/ [MeV] : Sp, = Sp^ = 2.6670 (2) 

Sf] = Ef] = 2.5710 

We use KB3pf for the case when the single-particle energies are set to their overall average: 

KB3p/ [MeV] : ep = ef = 2.6 (3) 

Due to the near degeneracy of the single-particle energies of the KB?>p_f interaction (0) , the results for the KB3pf 
case are very similar to those for KBip_f . 



III. COMPUTATIONAL PROCEDURES 

The computational procedures and tools used in the analysis are described in this section. In brief, the Hamiltonian 
and other matrices are calculated using an m-scheme shell model code |ll) while the eigenvectors and eigenvalues are 
obtained by means of the Lanczos algorithm |l2[. All the calculations are done in the full fp -shell model space. 

First, the Hamiltonian H for each interaction {KB3 (|l|), KB3p-f (g), and KBSpf (||)) is generated. Then the 
eigenvalues and eigenvectors are calculated and the yrast states identified. Next, the matrix for the second order 
Casimir operator of SU{3), namely C2 = ji^L"^ + Q ■ Q), is generated using the shell model code and a moments 
method |13| ] is used to diagonalize the C2 matrix by starting the Lanczos procedure with specific eigenvectors of H 
for which an SU{3) decomposition is desired. Finally, B{E2) values in e^frn^ units are calculated from one-body 
densities using Siegert's theorem with a typical value for the effective charge [|lj], qeff = 0.5, so Cp = {l + qeff)e = 1.5e 
and e„ — {qeff)e = 0.5e. 

Although the used procedure can generate the spectral decomposition of a state in terms of the eigenvectors of 
C2 of S'[/(3), this alone is not sufficient to uniquely determine all irreducible representation (irrep) labels A and /i of 
5(7(3). For example, C2 has the same eigenvalue for the (A,(u) and (/U, A) irreps. Nevertheless, since for the first few 
leading irreps (largest C2 values) the A and /i values can be uniquely determined [Esj this procedure suffices for our 
study. 

Usually, when considering full-space calculations, a balance between computer time and accuracy has to be consid- 
ered. While the Lanczos algorithm ||l2[ is known to yield a good approximation for the lowest or highest eigenvalues 
and eigenvectors, it normally does a relatively poor job for intermediate states. This means, for example, that higher 
states, in particular high total angular momentum states, may be poorly represented or, in a worst case scenario, not 
show up at all when these states are close to or beyond the truncation edge of the chosen submatrix. An obvious 
way to maintain a good approximation is to run the code for each Mj value, that is, Mj = 0, 2, 4, 6. . . . However, 
this might be a very time consuming process, but nonetheless one which could be reduced significantly if only a few 
Mj values are used for each run. For the calculations of this study, we used Mj = 0,6, 10, and 14. To maintain 
high confidence in the approximation of the intermediate states which have J = 2,4, 8, 12, ... we required that they 
be within the first half of all the states produced. The code was set up to output 29 states. A further verification on 
the accuracy of the procedure is whether the energies of the same state calculated using different Mj runs are close 
to one another. For example, as a consistency check the energy of the lowest J = 6 state in the Mj = run was 
compared to the energy of the same state obtained from the Mj = 6 run. 

IV. RESULTS 

Results for the SU{Z) content of yrast states and their B{E2) values for representative /p-shell nuclei are reported 
in this section. We focus on ^"'Ti, ^^Ti, ^^Ti, and ^^Cr because these are fp-sheW equivalents of ^"A^e, '^'^Ne, '^'^Ne, 
and '^'^Mg, respectively, which are known to be good SU{3) sd-shell nuclei. Furthermore, data on these nuclei are 
readily available from the National Nuclear Data Center (NNDC) llq| and full fp-sheW calculations are feasible [ pT[ . 
The model dimensionalities for full-space calculations increase very rapidly when approaching the mid-shell region; 
those for the cases considered here are given in Table I. 

Nucleus Mj = Mj = 6 Mj = 10 Mj = 14 

^Wi 1080 514 30 — 

^'^Ti 43630 32297 4693 134 

^^Ti 317972 278610 57876 3846 

*^Cr 492724 451857 104658 8997 

TABLE L Space dimensions for m-scheme calculations in full /p-shell model space. The computer code uses even parity and 
even isospin basis states with no restrictions on the total angular momentum J except for Alj — case where the computer 
code selects only states with even J values. 



In the following, we use four different graphic representations to illustrate our results. The first set, Figs.|| and ||, 
demonstrates the recovery of the SU{3) symmetry as the single-particle spin-orbit interaction is turned off, that is, 
in going from the KB3 to the KB3p-f interaction. Corresponding results for the KB3pf interaction are not given 
since they are similar to the KB3p-f results. In each graph, C2 values of SU{3) are given on the horizontal axis with 
the contribution of each SU{3) state on the vertical axis. The bars within each cluster are contributions to the yrast 
states starting with the ground state (J — 0) on the left. Hence the second bar in each cluster is for the J = 2 yrast 
state, etc. 
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Values of the second-order Casimir operator of SU(3) 

FIG. 1. Strength distribution of C2 of SU{3) in yrast states of '^'^Ti for realistic single particle energies with Kuo-Brown-3 
two body interaction {KB3). 
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FIG. 2. Strength distribution of C2 of 5f/(3) in yrast states of "^^Ti for degenerate single particle energies with Kuo-Brown-3 
two body interaction {K BSp^f) . 

We have chosen ^^Ti for an in-depth consideration of the fragmentation of the C2 strength in yrast states. The 
results on the nondegenerate KB3 interaction are shown in Fig.|l|. In this case the highest contribution (biggest 
bar) is more than 50% which corresponds to a C2 value of 114 for the J = 12 state. The C2 = 114 value is for 
(A, fi) — (8, 2) which is two SU{3) irreps down from the leading one, (A, /x) = (12, 0) with C2 = 180. The leading irrep 
only contributes about 10% to the J = 12 yrast state. The contribution of the next to the leading irrep, C2 = 144 for 
(A, fi) = (10, 1), is slightly less than 40%. Thus, for all practical purposes, the first three irreps determine the structure 
of the J = 12 yrast state. This illustrates that the high total angular momentum J states are composed of only the 
first few SU{S) irreps. This is easily understood because high J values require high orbital angular momentum L 
which are only present in SU{3) irreps with large C2 values. The high J states may therefore be considered to be states 
with good SU{3) symmetry. However, this is not the case with the ground state of '*'*Ti which has very important 
contributions from states with C2 values 60, 72, 90, 114, 144, and 180 with respective percentages, 7.5, 25, 10, 21, 8, 
and 21%. This shows that the leading irrep is not the biggest contributor to the J — ground state; there are two 
other contributors with about 20%, the third (C2 = 114) and seventh (C2 = 72) SU{3) irrep. 



When the spin-orbit interaction is turned off, which yields nearly degenerate single-particle energies since the single- 
particle orbit-orbit splitting is small, one has the KBip.f interaction and in this case the structure of the yrast states 
changes dramatically, as shown in Fig||. From Fig.g one can see that the leading irrep plays a dominant role as its 
contribution is now more then 50% of every yrast state. As in the previous case, the high total angular momentum J 
states have the biggest contributions from the leading irrep, for example, more than 97% for J — 12, 91% for J — 10, 
and 80% for J = 8. The ground state is composed of few irreps with C2 values 72, 114, and 180, but in this case 
the leading irrep with C2 = 180 makes up more than 52% of the total with the other two most important irreps 
contributing 21% [C2 = 72, (A,//) = (4,4)] and 23% [C2 = 114, (A, a*) = (8,2)]. 

An alternative way to show these results is given in FigJ| and Fig.H. These figures show the centroid, width, 
and skewness of the C2 distributions. The J values are plotted on the horizontal axis with the centroids given on 
the vertical axis. The width of the distribution is indicated by the length of the error bars which is just the rms 



(C2 — (^2)) ), from the average value of the second-order Casimir operator (C2). The third 



deviation, AC2 



central moment, SCo 



above, AC2 4- ^^, and below, AC2 — ^^, the average value. 



(C2 — {C2)) \ which measures the asymmetry, is indicated by the length of the error bar 
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FIG. 3. Average C2 values for KB?> and KB3p-f interactions in ^^Ti. 
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FIG. 4. Average C2 values for KB3 and KB3p-f interactions in ^^Ti. 



Note that the recovery of the leading irrep when the spin-orbit interaction is turned off is clearly signaled not only 
through an increase in the absolute values of the first centroid (C2) but also through the skewness SC2- For example, 
in *'^Ti with the KB3 interaction (spin-orbit interaction turned on) the ground state J = has (C2) = 110 and 
skewness SC2 = 33. This changes for the KB3p_f interaction to (C2) = 139 and a skewness of SC2 — ~37, as shown 
in Fig.H The equivalent of the '^'^Ti graph for the '^^Ti case is shown in Figjj. As for the ^'^Ti case, the results show 



the recovery of the SU{3) symmetry in ^^Ti when the single-particle spin-orbit interaction is turned off. 

We now turn to a discussion of the coherence nature of the yrast states. First notice that the widths of the 

distributions as defined by AC2 = \ ({C2 — (C'2)) ) are surprisingly unaffected (Fig.^ and Fig.^ by turning the 



spin-orbit interaction on and off. This effect occurs in all cases studied: ^^Ti, ^^Ti, "'^Ti, and '^^Cr. The more 
detailed graphs, Fig.Q and Fig.g, offer an explanation in terms of the fragmentation of the C2 distribution. As can 
be seen from these graphs, the irreps that are presented in the structure of a given yrast state in the presence of 
the spin-orbit interaction (Fig.nl) remain present, even though with reduced strength, in the structure of the state 

when the spin-orbit interaction is turned off (Fig.g). As a consequence, AC2 = \ ({C2 — (C'2)) ) which measures the 



overall spread of contributing irreps, is more or less independent of the spin-orbit interaction. One can see a sharp 
decrease in the width of the distribution only for high spin states like J = 12 in the graph for '^^Ti in Fig.^ 

The third type of graph, FigJq, demonstrates the coherent nature of the states within the yrast band. The three 
graphs shown give the spectrum of the second-order Casimir operator C2 of SU{i) for the J = 0, 2 and 4 yrast states 
in '^^Cr. The axes are labelled the same way as in Figs.n^andH, but in this case all bars are for a single yrast state. In 
this figure there are three peaks surrounded by smaller bars that yield a very similar enveloping shape for the given 
yrast states. The fragmentation and spread of C2 values is nearly identical for these states with no dominant irrep, 
indicative of severe S'f/(3) symmetry breaking. 
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FIG. 5. Coherent structure of the first three yrast states in "^^Cr calculated using realistic single particle energies with 
Kuo-Brown-3 two body interaction {KB3). On the horizontal axis is C2 of SU{3) with contribution of each 5(7(3) state to the 
corresponding yrast state on the vertical axis. 



Graphs for the KB3p-f case, when the spin-orbit interaction is turned off, are not shown since the results are 
similar to the results for ^^Ti shown in Fig.^. For example, when the spin-orbit interaction is on (KB3) the leading 
irrep for "^^Cr has a C2 value of 396 and this account for only around 10% of the total strength distribution (see 
Fig.|5|), but when the spin-orbit interaction is off (KB3p_/) the leading irrep is the dominant irrep with more than 
55% of the total strength. 

The last type of graph, Figs.g, 0, and g shows relative B{E2) values, that is, B{E2) strengths normalized to the 
B{E2 : 2+ — > 0+) value. For isoscalar transitions the relative B{E2) strengths are insensitive to the effective charges 
which may be used to bring the theoretical B{E2 : 2"*" —^ 0"*") numbers into agreement with the experimental values. 
Whenever an absolute B{E2 : 2+ -^ 0+) values are given they are in e^ frn^ units and the effective charges are 1.5e 
for protons and 0.5e for neutrons (<?e// = 0.5). 

The first graph on relative B{E2) values ( FigJGl) recaps our results for ^^Ti. Calculated relative B{E2) values for 
'^'^Ti corresponding to the spin-orbit interaction turned on (KB3) and spin-orbit interaction off (KB3p_/) are very 
close to the pure SU{3) limit. The agreement with experiment is very satisfactory except for the 4+ -^ 2+ and 
8+ -^ 6+ transitions. However, the experimental data ||l6[ on 8+ -^ 6+ transition gives only an upper limit of 0.5 
pico-seconds to the half-life. We have used the worse case, namely a half-life of 0.5 ps, as a smaller value would 
increase the relative B{E2) value. For example a half-life of 0.05 ps will agree well with the relative B{E2) value 
for the KB3p-f interaction. This example supports the adiabatic mixing which seems to be present for all the yrast 
states of ^-^Ti. 
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FIG. 6. Relative B{E2) values ( b(e2:2+^o+) ) ^^^ '^*'^^- ^he B{E2 : 2+ ^ 0+) transition values are 122.69eVm^ from 
experiment, 104.82e^/m'' for the KB3 interaction, and ISS.SSe^/m* for the KB3p-f case. 



Fig.0 shows B{E2) values for ^^Ti. In this case there are deviations from adiabatic mixing for the 6+ -^ 4+, 
10+ -^ 8+, and higher transitions. Two experimental data sets are shown in Figj^: data from the NNDC is denoted 
as Exp_(NNDC), and updated data on 2+ -^ 0+ and 4+ -^ 2+ transitions from pSJ is denoted as Exp_(Updated) . 
For ^^T« the agreement with the experiment is not as good as for ^'^Ti, however the experimental situation is also 
less certain. However, the coherent structure is well demonstrated for the first three yrast states 0+, 2+, and 4+ via 
relative B{E2) values for the KB3 and KB3p_f interactions which are very close to the SU (3) limit. 
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FIG. 7. Relative B{E2) values ( b{E2-2+Zo+) ) ^°^ ^^^^- ^he B{E2 : 2+ -^ 0+) transition values are 199.826^/771** from 
experimental data, 181.79e^/7n^ from updated experimental data, 208e^ fm'^ for KB3 interaction, and 299.83e^/7n'' for KB3p-f. 
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FIG. 8. Relative B[E2) values ( bIe2'2+Zo+] ) ^^ ''*^*- ^he B{E2 : 2+ -» 0+) transition values are 144. 23e^ fm'^ from ex- 
perimental data, 155.56^/777^ from updated experimental data, 202.46^/771* for KB3 interaction, and 445.326^/777* for KB3p-f. 

We conclude this section by showing the recovery of the SU{3) symmetry; this time via relative B{E2) values as 
shown for ''^Ti in Fig.g. In Fig.||we see that for the degenerate single particles case (KB3p_/) the first few transitions 
have relative B{E2) values which follow the SU{3) limit very closely. On other hand, the interaction involving spin- 
orbit splitting (KB3) is far from the SU{3) limit. The B{E2 : 4+ -^ 2+) transition is strongly enhanced due to the 
adiabatic mixing which is missing in the higher than J = 4 yrast states. 



V. CONCLUSION AND DISCUSSION 



The results reported in this paper show that the single-particle spin-orbit splitting is the primary interaction 
responsible for breaking of the 5^/(3) symmetry for nuclei in the lower /p-shell. When the spin-orbit splitting 
is reduced, as in the KBSp_f case, the importance of SU{3) as seen through the dominance of the leading irrep 
represented in each yrast state is revealed. It is important to note in this regard that the p- and /-shells are nearly 
degenerate, which implies a small P splitting. 

Although the SU{3) structure of the states is lost in the lower /p-shell, the results also show the mixing of SU{3) 
irreps that occurs displays enhanced B(E2) strengths. This adiabatic mixing results in a coherent structure that is 
represented in all yrast states for the '^^Ti case, while for the other nuclei studied this coherence breaks down after 
the first few yrast states. In particular, even though the yrast states are not dominated by a single SU{3) irrep, the 
B{E2 : 4+ -^ 2+) values remain strongly enhanced with values close (usually within 10-20%) to the SU{3) symmetry 



limit. 
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